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We show by extensive numerical simulations and analytical variational calculations that elongated 
binary non-miscible Bose-Einstein condensates subject to periodic modulations of the radial confine- 
ment exhibit a Faraday instability similar to that seen in one-component condensates. Considering 
the hyperfine states of ^^Rb condensates, we show that there are two experimentally relevant sta- 
tionary state configurations: the one in which the components form a dark-bright symbiotic pair 
(the ground state of the system) , and the one in which the components are segregated (first excited 
state). For each of these two configurations, we show numerically that far from resonances the 
Faraday waves excited in the two components are of similar periods, emerge simultaneously, and do 
not impact the dynamics of the bulk of the condensate. We derive analytically the period of the 
Faraday waves using a variational treatment of the coupled Gross-Pitaevskii equations combined 
with a Mathieu-type analysis for the selection mechanism of the excited waves. Finally, we show 
that for a modulation frequency close to twice that of the radial trapping, the emergent surface 
waves fade out in favor of a forceful collective mode that turns the two condensate components 
miscible. 
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I. INTRODUCTION 



The excitation of surface waves through parametric 
resonances is one of the oldest pattern-forming processes 
that goes back to Ernst Chladni's "beautiful series of 
forms assumed by sand, fillings, or other grains when 
lying upon vibrating plates" that "are so striking as to 
be recalled to the minds of those who have seen them 
by the slightest reference", Hans Christian 0rsted's ex- 
periments with lycopodium light powders, and Michael 
Faraday's "crispations" seen in "fluids in contact with vi- 
brating surfaces" [1]. The prototypical example of para- 
metric wave excitation is that of a shallow disc of a liquid 
rigidly oscillated in the vertical direction. In this setting 
the acceleration periodically modulates the effective grav- 
ity and for drives of sufficiently large amplitudes a surface 
wave instability occurs with frequency one half that of the 
drive Q. Such surface waves are termed Faraday waves, 
and the corresponding pattern-forming phenomenon has 
been seen in numerous Newtonian and non-Newtonian 
fluids, colloidal suspensions, ferromagnetic bodies, and, 
more recently, in superfluids. 

After a series of inceptive studies on extended paramet- 
ric resonances in confined superfluids [3-8], the experi- 
mental observation of Faraday waves in ^He cells [9| and 
^''Rb cigar-shaped Bose-Einstein condensates (BECs) 
[lo| catalyzed the interest in the nonlinear dynamics of 
parametrically-driven ultracold gases [ll|. In addition 
to the ever-present collective oscillation modes of BECs 
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[HI , their parametric driving can generate soundlike den- 
sity waves, which are analogous to Faraday surface waves. 
BEC systems also exhibit resonances, and when the driv- 
ing frequency is close to one of them, forceful resonant 
waves can develop, eventually taking over the dynamics 
of the system. Furthermore, inherent nonlinearity of such 
systems brings about a diverse set of related dynamical 
phenomena and leads to their complex interplay. 

Surveying the recent literature for bosonic systems 
one notices the theoretical investigations into the soli- 
ton management in periodic systems [l3|, fl3, t he emer- 
gence and suppression of Faraday patterns jl5Ml9| . the 
parametric excitation of resonances [20] and "scars" in 
BECs [21], spatially and teuiporally driven atomic in- 
teractions in optical lattices [22|, |23| , quantized vortices 
induced by spatio-temporally modulated interaction [2^ , 
stability and decay of Bloch oscillations in the presence 
of time-dependent nonlinearity [25|, [26|, and, quite in- 
terestingly, the removal of excitations in BECs subject 
to time-dependent periodic potentials [27]. On a related 
topic, the formation of density waves has been predicted 
for expanding condensates [29| and the spontaneous 
formation of density waves has been recently reported 
for antiferromagnetic BECs [30|. Faraday waves have 
been analyzed in detail also in superfiuid fermionic gases 
im and it has been shown that the collective modes 
of ID fermionic systems can be amplified by parametric 
resonances to the extent of observing a clear spin-charge 
separation [33]. 

Parametric resonances in BECs are usually achieved 
through periodic modulation of the frequency of the trap- 
ping potential, as is the case in Ref. [10], but the re- 
cent experiments on the collective modes of a trapped 
''Li BEC through periodic modulation of the scattering 
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length [34] have opened a new direction for both experi- 
mental and theoretical [35,] investigations. Simultaneous 
modulations of the strength of the confining potential and 
of the scattering length, in particular, are largely unex- 
plored. They can give rise to new recipes for pattern 
formation and, possibly, to new types of patterns. 

In this paper we study the excitation of waves through 
periodic modulations of the radial confinement of binary 
non-miscible BECs and show that these condensates ex- 
hibit a Faraday instability similar to that seen in one- 
component systems. Considering the hyperfine states 
of realistic ^^Rb condensates which can be readily pro- 
duced such as 1), |2,0) and 1), |2, 1) pairs, 
we show that there are two distinct experimentally rel- 
evant stationary configurations: one where the compo- 
nents form a dark-bright symbiotic pair, which is the 
ground state of the system, and one where the compo- 
nents are segregated, which is the first excited state of 
the system. Far from resonances we show numerically 
for each configuration that the excited waves are of sim- 
ilar periods and emerge simultaneously, and analytically 
find the dispersion relation using a variational treatment 
in conjunction with a Mat hieu- type analysis. Finally, the 
resonant excitation of collective modes is analyzed in de- 
tail. 

The rest of the paper is structured as follows. In 
Sec. ini we describe the numerical treatment of the cou- 
pled Gross-Pitaevskii equations (GPEs) that describe the 
T = dynamics of the condensate and introduce the two 
types of stationary configurations, while in Sec. lIIII we de- 
rive the corresponding variational equations and the as- 
sociated dispersion relations. In Sec. [IVlwe present our 
numerical and analytical results for realistic ^^Rb con- 
densates, along with suggestions for future experiments. 
The last section contains conclusions and outlook. 



II. STATIONARY STATES 

The dynamics of binary condensates at T = is 
governed by the time-dependent version of the coupled 
GPEs, which read 

^Ns-jU\^3-j\'^j. (1) 

where j G {1,2}, Nj is the fixed number of atoms in 
each component, Uj = 47raj, U = 47ra, with aj being the 
intra-component scattering lengths, while a is the inter- 
component scattering length. Here we use numerical val- 
ues similar to the experimental ones from Refs. [38-40], 

ai = 100.4 ao, a2 = 98.98 ao, a = ai , (2) 

where ao is the Bohr radius. Small variations in scat- 
tering lengths yield similar results, but we have cho- 
sen these particular values since they correspond to a 
clearly non-miscible configuration at the mean-field level. 



thereby emphasizing the forcefulness of the miscibility 
transition, which represents one of our main motivations, 
as we will see in Sec. |lVl For simplicity, we use natural 
units h = m = 1 throughout the paper, and component 
normalization 

j dv\i^i{v,tt = 1. (3) 

The stationary states of binary Bose-Einstein conden- 
sates are obtained by solving the time-independent ver- 
sion of Eqs. ([1]), in which terms with time derivatives 
of wave functions are replaced by the terms containing 
chemical potentials fij of components, 

+7V3_^-/7|^3-,N, = ^j^j- (4) 

In this paper we consider the two hyperfine states of ^'^Rb 
(hereafter referred to as states A and B) in an external 
potential of the form 

V{v,t) = \nl{t)p' + \^llz\ (5) 

where = such that the system is cylindrically- 

symmetric, i.e. tljj{r^t) = ipj{p^z^t)^ and GPEs are effec- 
tively two-dimensional. We also assume that the system 
is highly elongated and strongly confined in the radial 
direction, ftp{t) ^ Qz- The above time-dependent sys- 
tem of GPEs as well as its time-independent coun- 
terpart Eq. (jl]), can be solved numerically using various 
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FIG. 1: (Color online) Typical imaginary-time propagation of 
the total energy for a two-component condensate system with 
A^i = 2.5 • 10^ atoms in the state A and iV2 = 1.25 • 10^ in 
the state B, Qp = 160 x 27r Hz, = 7 x 27v Hz. The dashed 
black curve corresponds to the evolution from initial wave 
functions set to two identical Gaussians, which immediately 
yields the ground state, while the full red curve corresponds 
to evolution from the initial wave functions set to two well- 
separated Gaussians, yielding the first excited state of the 
system, which eventually decays to the ground state. 
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approaches [4l|-l45|. Here we choose the efficient spht- 
step Crank-Nicolson approach developed by Adhikari and 
Muruganandam in Ref. [46|. We have implemented this 
algorithm in the C programming language, and use it for 
all numerical simulations presented here. 

Binary condensates exhibit a wide range of interest- 
ing and experimentally relevant configurations, which 
go from the one-dimensional soliton pairs tabulated in 
Ref. [43jto the more exotic two-dimensional vortex- 
lattices jlil, and vortex-bright-soliton structures [49|. 
From the mean-field theory it follows that binary con- 
densates are miscible if the condition U < ^/UiU2 is sat- 
isfied [50]. For the values from Eq. ([2]) this condition is 
not satisfied, and the system is non- miscible (ssj . 

The stationary states are computed numerically by 
imaginary-time propagation until we achieve the conver- 
gence of wave functions and physical quantities of the sys- 
tem, chemical potentials and mean-square radii of com- 
ponents, and the total system energy. Fig. [T] illustrates 
the results obtained using this approach for different ini- 
tial conditions. It shows total energy of the system as 
a function of the imaginary time. Using the various ini- 
tial conditions, we were able to numerically calculate two 
relevant stationary configurations: the one in which the 
components form a dark-bright symbiotic pair (ground 
state), and the one in which the components are spa- 
tially well segregated. From Fig. [T] we can see that the 
segregated state is a first excited state of the system, 
as it decays to the ground state after sufficiently long 
imaginary-time propagation. 

The first type of stationary solutions supported by 
the GPEs, the ground state, consists of a symmetric, 
Thomas- Fermi-type density profile with a hole in the cen- 
ter in one component, and a well- localized density peak 
positioned in the center of the trap in the other compo- 
nent. Fig. [2] shows typical longitudinal density profiles 



nj{z) = dp 27rp 
Jo 



(6) 



of the two components for the ground state, obtained 
through the imaginary-time propagation, starting from 
the two identical Gaussian initial states. 

We note the apparent similarity of the obtained dark- 
bright symbiotic solutions and the dark-bright soliton 
molecules, known to exist in homogenous systems. These 
soliton molecules were originally predicted in a nonlinear 
optics setting [si], [52l and were first observed in pho- 
torefractive crystals [53| . Following their exposure to the 
BEC community [si] , they were observed experimentally 
using the ^'^Rb condensates [55|, and were subsequently 
addressed theoretically in Refs. [39|, HH, [H^- In homo- 
geneous condensates the dark component effectively acts 
as a trapping potential for the bright component through 
the nonlinear interaction, and this mechanism is well pre- 
served in inhomogeneous systems, albeit the dark compo- 
nent is cut-off away from the center of the trap and the 
bright component tends to be more narrow due to the 
additional confinement by the trap. However, symbiotic 
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FIG. 2: (Color online) Typical longitudinal density profiles 
for a condensate in the ground state (dark-bright symbiotic 
pair). The parameters are the same as in Fig. [U the full red 
curve shows the density profile for atoms in the state A, the 
dashed black line for atoms in the state B. 
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FIG. 3: (Color online) Typical longitudinal density profiles 
for a condensate in the first excited state (segregated state). 
The parameters are the same as in Fig. [U the full red curve 
shows the density profile for atoms in the state A, the dashed 
black line for atoms in the state B. 



solutions from Fig. [2] do not represent inhomogeneous 
counterparts of soliton molecules, since they are purely 
real, and the first component does not exhibit a jump of 
TT in the phase between its left and right part. 

The second type of stationary solutions supported by 
the GPEs, the first excited state, is presented in Fig. [3l 
It consists of two opposing, asymmetric Thomas-Fermi- 
type density profiles, with the negligible overlap. This 
segregation of the components stems from the repulsive 
character of the interaction and the non-miscible nature 
of the considered binary condensate. Naturally, there ex- 
ist a large number of further excited states of the system 
that can be calculated numerically as well, but we focus 
only on the two relevant states, shown in Figs. [2] and [3l 

Experimentally, such states can readily be realized us- 
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ing current technology. While our analysis considers 
elongated traps approaching a one-dimensional limit, we 
note that the state in Fig. [2] has already been observed 
in three-dimensional traps where a magnetically trapped 
BEC of ^^Rb atoms was prepared in a mixture of the 
1 1,-1) and |2,1) hyperfine states and was observed to 
assume a ball-shell structure [58', '59^ . By introducing dif- 
ferential trap shifts for the two components as in Ref. [SSj , 
the segregated state in Fig. [3] can readily be generated as 
well. Alternatively, optical traps can be employed to con- 
fine atoms in a state independent way, and an additional 
magnetic gradient can be used to separate components of 
a mixture. Using a miscible mixture, this phase separa- 
tion in an optical trap has been generated e.g. in Ref. 1401 
and can also be extended to optical lattice systems [6^. 
Using a microwave sweep in a high-bias magnetic field, 
the miscible mixture in the latter experiments can be 
converted into an immiscible mixture by changing the 
hyperfine state of one of the components. 

At the end, let us also note that the stationary char- 
acter of the two numerically calculated states presented 
here was cross-checked by performing a real-time propa- 
gation for long time intervals (much longer than those 
considered in the rest of the paper). While in the 
imaginary-time propagation the first excited state even- 
tually decays to the ground state, when real-time prop- 
agation is performed both stationary states are found to 
be stable for all practical purposes. 



III. VARIATIONAL TREATMENT 

In this section we develop variational analytic ap- 
proach suitable for study of the emergence and char- 
acterization of Faraday patterns in non-miscible binary 
condensates, induced by harmonic modulation of the ra- 
dial trapping frequency. For each of the two stationary 
states calculated in Sec. |TT1 we propose a suitable vari- 
ational ansatz for component wave functions, approxi- 
mately solve the ensuing equations, and analytically de- 
rive expressions for periods of Faraday waves induced 
far from resonances. The obtained analytical results are 
compared with the numerical results in Sec. IIVI 



A. Symbiotic pair state 

To construct a suitable variational ansatz for compo- 
nent wave functions in the case of the symbiotic pair 
state, we consider an equivalent one-component scenario, 
in which the bright state ip2 evolves in the combined 
effective field created by the dark state and by the 
trapping potential. To further simplify the problem, in 
this subsection we will also assume that the values of 
three scattering lengths from Eq. (j2j) are equal, namely 
ai= a2 = CL = a = 100.4 ao, or Ui = U2 = U = /7, which 
is justified since the values measured experimentally are 
quite close [ssj . 



In this setting, the effective one-component Lagrangian 
density for the bright state ?/^2 is given by 



dt 



■ d^p2 
dt 



^V{p,Z,t)\l^2 



UN2 I , ,4 



-^UN^\^,f\^2f. 



(7) 



where we have chosen the following ansatze for compo- 
nent wave functions: 



ilji{p,z,t) = A/^exp 



-ip^a^(t) 



1 — exp 



2wl) 



z 



2wl{t) 

X [1 + {u{t) + iv{t)) cos kz\ 



(8) 



(9) 



with straightforward interpretation for the variational 
parameters: the radial and the longitudinal bright state 
widths Wpit) and the phase o?{t)^ and the complex 
amplitude u{t) + iv{t) of the Faraday wave in the bright 
component, with the period 27r//c. Note that, in order 
to keep the analytics reasonably simple and tractable, 
we assume that the longitudinal condensate width Wz is 
constant, and include the surface wave only in the bright 
component, in a manner similar as in Ref. j6ll-[63j. As 
the role of the dark component is mainly that of an addi- 
tional trapping potential, these simplification of ansatze 
have little impact on the final results. The normalization 
factors Mj are calculated from normalization conditions 



dz / dp 27rp|V^i(p,z,t)r = 1, (10) 



dp2np\^2{p.z,t)\^ = 1, (11) 



/OO nC 
dz / 
-00 ^0 



where 2L is the longitudinal spatial extent of the dark 
state component, obtained by solving the stationary sys- 
tem of GPEs and assumed to be constant. 

After inserting the ansatze (|8|) and (|9]) into the La- 
grangian density ([7j), we calculate the Lagrangian and 
derive the following variational equations for the param- 
eter functions Wp(t)^ <^(^), u{t) and v{t)'. 



2wpa^ 



2 

-2a^ 

2 



+ 



U {3VSLN2^gaV^Wz) 



u = 



2471^ w'^Wz [2L- {VS-l) 

Pu UN2U 

2 



^/2ii'2WzW 



(12) 

(13) 

(14) 
(15) 
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where = {12^6V2-SV6)Ni-3{A- V2)N2. In addi- 
tion to these ordinary differential equations, we have an 
algebraic equation for namely 



= 1 



V2 



UwzN27r 2 



1 



L2 



22 



{2V2-l)V^Wz 



(16 



L 

2i)y/l\Wz 



(16) 



irgzW^ 



where gz = 2(9^2 - 16^3 + 4^6 + 6) A^i - 3(9^2 - 8)7V2 
and Wpo = Wp{0). The tacit assumption that the conden- 
sate is of constant longitudinal extent 2L will be further 
justified in the next section, where we present the full 
numerical results. 

Let us first emphasize that, as will be shown later nu- 
merically, the dynamics of the bulk of the condensate, in 
the first approximation, is not impacted by that of the 
surface wave. Eqs. (p!2l) and (p!3|) are, in fact, similar to 
those derived in Ref. [64] for the collective dynamics of 
low-density one-component condensates, while Eqs. ([H]) 
and (p!5|) resemble those derived in Refs. [i5|,[l7] for Fara- 
day waves in low- and high-density one-component con- 
densates. Second, we stress that for a modulated radial 
trapping ftp{t) = ftpo • (1 + esincjt), Eqs. ([12]) and ([13]) 
exhibit a series of parametric resonances for cu = Qpo 
(self-resonance) and uj = 2Q.pQ/n^^ where n is an inte- 
ger Q. The widest resonance is that at = 2VtpQ (for 
n = 1), followed by the self-resonance at = Vtp^^ which 
has been evidenced in Ref. [io| for one-component ^^Rb 
condensates. The other (n > 2) resonances are very nar- 
row and are of little experimental interest, since for such 
low frequencies the excited waves have periods compara- 
ble to the longitudinal extent of the condensate and are 
therefore difficult to observe. 

Finally, let us note that, unlike in one-component con- 
densates, where Faraday waves emerge rapidly enough 
as to hide the resonant behavior at = 2VtpQ [lo|, in 
two-component systems the forceful resonant behavior is 
dominant, as we will show numerically in Sec. IIVI 

Far from resonances we can approximately calculate 
the radial width as 



Wp ^ Vtp 



U{3VSLN2^gaV^Wz) 
Uiriwz [2L - (VS-l) v/tF^^] 



(17) 

which stems from Eqs. ([T2)) and (fT3|) by neglecting Wp(t)^ 
and casts the equation for u{t) into the form 

u{t) + u{t) [a{k, uj) + eh {k, uj) sin 2r] = , (18) 

A dimensionless time r is introduced as ut = 2r, and 



^4 ^2 

a{k,uj) — + — Asym, 



6(/c, uo) 



A, 



sym 



(19) 
(20) 
= ,(21) 



Qah 



= 2L- (2\/2 - l)V7r^ 



TTWz 



(22) 



For small positive values of the radial modulation am- 
plitude e and positive values of the function b{k^uj), 
Eq. ([T8|) has solutions of the form exp(±z/ir) sin ^/ar 
and exp(±i/ir) cos ^/ar, where lm[/j] consists of a series 
of lobes positioned around the solution of the equation 
a{k^uj) = n^, with n being an integer [65|. The lobe cen- 
tered around a{k^uo) = 1 is the largest, and it yields the 
most unstable solutions [HI, [I^ , determined by 



As 



sym 



A2 



(23) 



As these solutions have a frequency of cj/2, half that 
of the parametric drive, they are usually referred to as 
Faraday waves in honor of Faraday's classic study [1]. 

Close to a resonance [66], the approximation used 
above for Wp{t) breaks down, and one cannot generally 
construct the explicit equation for u{t). We know, how- 
ever, that the instabilities appear due to the resonant 
energy transfer between the radial mode and the surface 
wave, which entails that Wp{t) and u(t) are of equal fre- 
quency. For small values of the modulation amplitude e 
this requires that the condition a{k^uj) = 2^ is satisfied. 



B. Segregated state 

To consider the case of a segregated excited state, we 
build the variational equations starting from the habitual 
Gross-Pitaevskii Lagrangian density [37, 64], written for 
a two-component condensate in the form: 



dt 



^UNiN2\^ir\^ 



(24) 



To variationally describe the wave functions of two BEC 
components, we use cylindrically-symmetric ansatze tai- 
lored around the usual Gaussian envelopes [64] that de- 
scribe the bulk of the condensate, to which we graft a 
surface wave l6ll-l63l. 



2wj{t) 
{u{t) +iv{t)) cos kz] e (( 



(25) 



where represents Heaviside step function, and normal- 
ization factors are determined from 



/oo rc 
dz 
-oo Jo 



dp 27r p\iljj{p,z,t)\ 



1 



(26) 



As far as the longitudinal components are concerned, 
the trial wave functions consist of two opposing half- 
Gaussians of equal widths and amplitudes, positioned in 
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the center of the trap. As there is no overlap between 
the two components and the period of the excited wave 
is smaher than the longitudinal extent of the condensate, 
for small- amplitude waves the Euler-Lagrange equations 
can be written as 



2wpa , 



2wl 
2~ 



2 



gu 



(27) 
(28) 

(29) 
(30) 



where g = NfUi + N^U2 and TV = TVi + 7V2. As in the 
previous subsection, in writing the equations above we 
have assumed that the condensate has a frozen longitu- 
dinal dynamics (apart from the dynamics of the grafted 
wave), and that the corresponding width Wz is constant. 
Its value is determined from the algebraic equation 



1 

2^ 



2 



g 



= 0. 



(31) 



where Wpo = Wp{0). As in the case of symbiotic states, we 
will consider modulation of the radial trapping frequency 
of the form Qp{t) = Qpo- {l-\-e sin ujt). Following a similar 
reasoning, we find that the system again exhibits a series 
of parametric resonances: a self-resonance at = l^po, 
and series of resonances for u = 2VtpQ/n^^ where n is 
an integer. As before, the strongest resonance is that at 
u) 21^po- 

Far from resonances, we can approximate the radial 
width as 



1 



1/4 



(32) 



while Eqs. (j29j) and ([3Q|) can then be conveniently recast 
in the form of a Mathieu equation. 



u{t) + u{t) [a(/c, uo) + eh (/c, uo) sin 2r] = 
where the time r is introduced as uot = 2r, 

a{k,uj) = — + ^Aseg, 



0, 



uo 



and 



A. 



2 ^seg , 



(33) 

(34) 
(35) 

(36) 



As before, the most unstable solutions correspond to 
a{k^uo) = 1, and are determined by 



A, 



^F,seg 



■seg 



A2 

seg 



(37) 



Close to resonances, following a similar reasoning as in 
the case of symbiotic pair solution, we conclude that the 
radial mode Wr{t) and the surface wave u{t) are of equal 
frequency, which, for small values of the modulation am- 
plitude e, is determined by solving a{k^uj) = 2^. 



IV. RESULTS AND DISCUSSION 

To investigate the emergence of Faraday waves, we 
solve numerically the coupled set of time-dependent 
GPEs ([T]) for the exact values of scattering lengths from 
Eq.(j2j) and study the dynamics of the longitudinal den- 
sity profiles of BEC components. 



nj{z,t) 



Jo 



dp 27ip\il)j{p,z,t)Y 



(38) 



as well as their Fourier spectra. Our main result is that, 
far from resonances, Faraday waves of similar periods ap- 
pear simultaneously in both BEC components for both 
considered initial configurations (symbiotic pair state, 
segregated excited state), and that these waves have al- 
most no effect on the dynamics of the bulk of the con- 
densate (surface waves represent only a small perturba- 
tion of the stationary state). For the self-resonance at 
uo = VtpQ^ we show that surface waves appear consider- 
ably faster than the Faraday ones, and that, similar to 
the one-component case reported in Ref. [10], such reso- 
nant waves have smaller period than the Faraday waves. 
To understand this latter feature, we recall from the pre- 
vious section that the instability now sets in due to the 
resonant energy transfer between the radial mode and the 
surface wave, which changes the frequency of the surface 
wave and consequently the observed period. Finally, we 
numerically study the forceful resonant dynamics of the 
system at = 2VtpQ^ which has not been seen previously 
in one-component condensates [lo| . This strong resonant 
behavior is found to facilitate the dynamical transition 
of the system from non-miscible to the miscible state. 



A. Symbiotic pair states 

In Fig. m we show the emergence of Faraday waves in 
real-time dynamics of a two-component condensate with 
Ni = 2.5 • 10^ atoms in the state A and N2 = 1.25 • 
10^ atoms in the state B, starting from a symbiotic pair 
ground state configuration. The magnetic trap has the 
parameters {Qpo^Qz} = {160 x 27r Hz, 7 x 27r Hz}, and 
we consider the modulation frequency uo = 250 x 27r Hz 
to be far from resonances. The Faraday waves are visible 
after 150 ms, and we show in Fig.[5]the Fourier spectrum 
of density profiles of the condensate at t = 200 ms. The 
spectrum exhibits several peaks, and the first two peaks, 
ki and /c2, common for both components, are related to 
the geometry of the system. The first peak corresponds 
to the period £1 = 27r//ci = 81.5 /im, the longitudinal 
extent of the system, while the second one corresponds 
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FIG. 4: Emergence of Faraday waves in a two- component 
BEG system in the real-time evolution of longitudinal density 
profiles: (a) ni(z, t) and (b) 712(^5 1). The system is initially in 
the symbiotic pair state, and is modulated with the amplitude 
e = 0.1 and frequency uj — 250 x 27r Hz. The system contains 
iVi = 2.5 • 10^ atoms in the state A and N2 = 1.25 • 10^ atoms 
in the state B, confined by the trap with r^^o = 160 x 27r Hz 
and = 7 X 27r Hz. 
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FIG. 5: (Golor online) Fast Fourier transform of density pro- 
files from Fig.[4]for two condensate egmpoppnts at t = 200 ms. 



to the period ^2 = 27r//c2 = 21.7 /im, the extent of the 
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FIG. 6: Emergence of resonant waves in a two- component 
BEG system in the real-time evolution of longitudinal density 
profiles: (a) ni(z,t) and (b) 712 (^,t). The system is initially 
in the symbiotic pair ground state, and is modulated with 
the self-resonant frequency uj — Qpo = 160 x 27r Hz; other 
parameters as in Fig. U] 



central dip in the fist component (or, equivalent ly, the 
extent of the second component). The periods of Faraday 
waves are determined by the peaks /ca^j, and have very 
close values, ^3^1 = 13.0 /im and ^3^2 = 12.5 /im. The 
dispersion relation (|23|) derived in Sec. IIII Al indicates a 
period of 12.0, /im, which is in excellent agreement with 
the numerical results. 

This demonstrates that the variational ansatze from 
Sec. IIII Al were well crafted. The good agreement is also 
partly due to our normalization of the dark component, 
where we have used the longitudinal extent of the conden- 
sate L = ^1 , obtained from the stationary solution of the 
full set of GPEs. Please note that the emergence of the 
Faraday waves does not impact the bulk of the conden- 
sate, and that its longitudinal extent is constant, which 
fully justifies our assumption in the variational model. 

In Fig. [6] we show the resonant dynamics of the con- 
densate for a self-resonance, uj = l^po = 160 x 27r Hz. 
In this case our numerical simulations indicate a narrow 
resonance where the excited surface waves do not im- 
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FIG. 7: Emergence of resonant waves in a two- component BEC system in the real-time evolution of longitudinal density profiles: 
ni{z,t) (left column) and 712(2;, t) (right column). The system is initially in the symbiotic pair ground state, and is modulated 
with the frequencies: cj = 300 x 2n Hz in (a) and (b); cj = 2Qpo = 320 x 27r Hz in (c) and (d); cj = 340 x 27r Hz in (e) and (f). 
Other parameters are the same as in Fig. [U 



pact the dynamics of the bulk of the condensate. The 
resonant waves develop faster than the Faraday waves, 
and already after 70 ms are clearly visible. However, in 
Fig. [3 we see that for a second resonance at cj ^ 2r^po, 
the excited surface waves appear even much earlier, after 
only 25 ms, but are quite short-lived, and the collective 
dynamics of the two components then takes over. This 
figure also illustrates that the second resonance is very 
broad, and covers the interval wider than 40 x 27r Hz, 



centered at 2Qpo — 320 x 27r Hz. 

We note, in particular, that precisely at the second 
resonance the two components explode into one another 
(due to the resonant energy transfer), thereby becoming 
effectively miscible. This resonant transition to miscibil- 
ity is specific to two-component systems, and has also 
been recently reported in binary dipolar BECs [g^ • Our 
analytical treatment of the surface waves indicates a pe- 
riod of 9.3 /im for the self-resonance, while the full GPEs 
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FIG. 8: Emergence of Faraday (a, b) and resonant (c, d) waves in a two- component BEG system in the real-time evolution of 
longitudinal density profiles: ni(z,t) (left column) and 712(^5 1) (right column). The system is initially in the segregated excited 
state, and is modulated with the amplitude e = 0.1 and frequencies: uj — 250 x 2n Hz in (a) and (b); uj — Qpo = 160 x 27r Hz in 
(c) and (d). The system contains Ni — 2.5 • 10^ atoms in the state A and A^2 = 1-25 • 10^ atoms in the state B, and is confined 
by the trap with Qpo = 160 x 27r Hz and = 7 x 27r Hz. 



numerically yield periods of 9.3 /im for the first com- 
ponent and 9.0 /im for the second component, which is 
again an excellent agreement. For the broad resonance at 
UJ ~ 21]p05 the variational treatment gives 4.9 /im, while 
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FIG. 9: (Golor online) Fast Fourier transform of density pro- 
files from Fig.[8]for two condensate components at t = 200 ms. 



numerically we find periods of 9.6 /im for both compo- 
nents. Here the agreement is only qualitative, due to a 
violent dynamics observed numerically, which cannot be 
fully captured by simple ansatze used in Sec. lEI A[ 

As symbiotic states are routinely obtained in ^'^Rb con- 
densates [39I, [Hll , the observation of the Faraday waves is 
definitely within the current experimental capabilities, as 
is the observation of self-resonant waves at = Qpo . The 
observation of resonant waves for uj close to the second 
resonance 2Qpo seems, however, unlikely, as these waves 
quickly fade out in favor of a forceful resonant dynamics 
that takes the condensate into the miscible regime and 
can, for longer timescales, turn the condensate unstable. 



B. Segregated states 

In Figs. [8^1 andlHb and we show the formation of Fara- 
day waves at = 250 x 27r Hz for a segregated initial 
state using the same experimental setup as in the previ- 
ous case. The waves again appear after roughly 150 ms 
and we show in Fig. [9] the Fourier spectrum of the longi- 
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FIG. 10: Emergence of resonant waves in a two- component BEG system in the real-time evolution of longitudinal density 
profiles: ni{z,t) (left column) and 712(^5 1) (right column). The system is initially in the segregated excited state, and is 
modulated with the frequencies: uj = 300 x 27v Hz in (a) and (b); a; = 2Qpo = 320 x 27r Hz in (c) and (d); a; = 340 x 27v Hz in 
(e) and (f). Other parameters are the same as in Fig. [51 



tudinal density profile at t = 200 ms. The spectrum's two 
peaks ki^i and ki^2 correspond to longitudinal extents of 
components, 81.5 /im and 54.3 /im, while the peak k2 cor- 
responds to the width of the interface area between the 
components of 21.7 /im. 

Periods of the Faraday waves are determined by the 
peaks ks^i and ks,2^ and have the values 11.6 jam for the 
first and 13.0 /im for the second component. The dis- 



persion relation derived in Sec. IIII Bl indicates a period 
of 16.8 /im, which overestimates the previous result by 
roughly 30% due to several reasons: the ansatze we use 
in Sec. IIIIBI for two components are symmetric around 
z = 0, they have the same normalization, and also they 
exhibit the Gaussian exponential cut-off away from the 
center of the trap. From Fig. [3] we already see that 
these are simplifications made only to ensure analytical 
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tractability of the variational calculation. We also stress 
that both components are in the Thomas- Fermi regime, 
hence their wave functions are well-localized and show 
a distinct quasi-singular cut-off at borders of the cloud, 
not the Gaussian-type tail. However, even with these 
simplifications, the variational approach is able to give 
reasonable estimate for the period of Faraday waves. 

We also numerically study resonant dynamics of the 
two-component BEG system for the self-resonance at 
uj = QpQ. In Figs. [Ht and [HJi we show real-time evolu- 
tion of density profiles for both components. The reso- 
nance is relatively narrow and the excited surface wave 
do not impact the overall dynamics of the condensate. 
The resonant surface waves appear after around 80 ms, 
and have the period of 11.2 /im in the first and 12.5 /im 
in the second component, while the variational approach, 
Eq. (|37|) . gives the period of 13.1 /im. The agreement is 
quite good, but only coincidental, since we know that the 
ansatze used in Sec. IIII Bl are oversimplified, and that we 
cannot expect more than an order of magnitude agree- 
ment. 

In the case of a second resonance at = 2Qpo^ the 
excited surface waves fade out in favor of the forceful 
resonant dynamics that turns the two components mis- 
cible. The resonant waves appear already after around 
25 ms, much faster than for the self-resonance. As in the 
case of the symbiotic pair state, the second resonance is 
much broader than the self-resonance, with the width of 
more than 40 x 27r Hz, as illustrated in Fig. [TOl Our ana- 
lytical treatment of the surface waves indicates a period 
of 6.6 /im, while the Fourier analysis of the full numeri- 
cal solution of GPEs yield periods of 10.2 /im in the first 
and 9.6 jam in the second component. As expected, the 
variational treatment gives a fair estimate of the period 
of resonant surface waves. 



V. CONCLUSIONS 

We have shown the emergence of Faraday waves in bi- 
nary non-miscible BEGs subject to periodic modulations 
of the transverse confinement. Considering the |1,— 1) 
and 1 2, 1) hyperfine states of realistic ^'^Rb condensates, 
we have shown that there are two types of experimen- 
tally relevant stationary configurations: a symbiotic pair 
ground state, where one component is trapped by the 
other, and a configuration where the components are well 
segregated, which represents a first excited state of the 
system. For both types of configurations we have shown 
by extensive numerical study that the surface waves in 
the two components emerge simultaneously, are of simi- 
lar periods and do not impact the dynamics of the bulk 
of condensate far from resonances. We have derived ana- 
lytically periods of the excited surface waves using first a 
variational treatment of the full system of coupled GPEs 
that reduced the dynamics of the condensate to a set of 
coupled ordinary differential equations, and then using a 
Mathieu-type analysis of these equations we have deter- 



mined the most unstable solutions. The obtained ana- 
lytical results for periods of surface waves are shown to 
give excellent estimates by comparison with the Fourier 
analysis of numerical solutions of the full set of GPEs. 
The resonant dynamics seen for modulations equal to 
the radial trapping frequency (self-resonance), as well as 
to twice this value (second resonance) were investigated 
numerically, and we have found that the emergent res- 
onant surface waves appear much faster than the Fara- 
day waves. Using the variational approach developed in 
Sec. Illli the periods of resonant surface waves can be 
also estimated, and are found to be in reasonable agree- 
ment with the numerical results. In the case of a second 
resonance, the emergent surface waves for both types of 
initial configurations are found to fade out in favor of a 
forceful resonant dynamics that turns the binary conden- 
sate miscible. This is an interesting phenomenon, and we 
believe that GPE framework captures its dynamics quite 
well, similarly to e.g. Ref. [68], where the condensate was 
split into two separated regions and then the collisional 
dynamics of the two halves was observed and found to 
be in excellent agreement with GPE predictions. 

Extending the present investigation to experimental 
setups with anisotropic transverse confinement (such as 
the one in Ref. presents one of interesting venues 

for future studies. As fully three-dimensional simulations 
are extremely time-consuming, a natural first step would 
be to extend the existing (cylindric ally symmetric) non- 
polynomial Schrodinger equations [69|, ItO] to cover the 
dynamics of binary BEGs subject to anisotropic trans- 
verse confinement. Also in this direction the study of pre- 
dominantly two-dimensional (pancake-shaped) conden- 
sates has the potential to uncover a rich variety of pat- 
terns, and one could even expect the formation of spatio- 
temporal chaos in such systems. Finally, another natural 
extension of this work would be to consider parametric 
resonances and the ensuing Faraday waves in miscible 
condensates. In this setting the dynamics of the waves 
in the two components can be mapped variationally onto 
a set of coupled Mathieu equations whose most unstable 
solutions are, however, not known analytically. 
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